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ABSTRACT 



Aims. We have searched for emission from the 557 GHz ortho-water line in the interstellar medium of six nearby starburst galaxies. 
Methods. We used the Odin satellite to observe the 1 10 - loi transition of o-H 2 in the galaxies NGC253, IC342, M82, NGC4258, 
CenA, and M51. None of the galaxies in our sample was detected. 

Results. We derive three sigma upper limits to the H 2 abundance relative to H 2 ranging from 2 x 10~ 9 to 1 x 10 -8 . The best of these 
upper limits are comparable to the measured abundance of H 2 in the Galactic star forming region W3. However, if only 10% of the 
molecular gas is in very dense cores, then the water abundance limits in the cores themselves would be larger by a factor of 10 i.e. 
2 x 10~ 8 to 1 x 10~ 7 . 

Conclusions. These observations suggest that detections of H 2 emission in galaxies with the upcoming Herschel Space Observatory 
are likely to require on-source integration times of an hour or more except in the very brightest extragalactic targets such as M82 and 
NGC253. 
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1. Introduction 



Recent space-based observations by S WAS (Melnick et al. 2000 ) 
and Odin (Frisk et al. |2003l l have revealed a low gas-phase abun- 
dance of water under most conditions in the interstellar medium 
(e.g. Snell et al. I2000al Olofsson et al. [2003T >. This low abun- 
dance is most likely due to the water molecules freezing out to 
form water ice on dust grains (e.g. Bergin et al. 2000). However, 
gas phase water abundances 100 to 10,000 times larger than the 
typical ambient cloud value have been seen toward hot cores 
with ISO (e.g. van Dishoeck & Helmich 1 19961 and in shock- 
heated outflows (Bergin et al. 120031) and supernova remnants 
(Snell et al. 120051) . The significant dependence of the gas-phase 
water abundance on environment makes it interesting to search 
for water emission lines in even more extreme environments, 
namely the nuclei of nearby starburst galaxies. 



Send offprint requests to: C. Wilson 

* Based on observations with Odin, a Swedish-led satellite project 
funded jointly by the Swedish National Space Board (SNSB), the 
Canadian Space Agency (CSA), the National Technology Agency of 
Finland (Tekes) and Centre National d'Etude Spatiale (CNES). The 
Swedish Space Corporation has been the industrial prime contractor and 
also is operating the satellite. 



In this paper, we present sensitive upper limits to the 1 io— lot 
emission line of ortho-HiO at 556.936 GHz in six nearby star- 
burst galaxies. The observations and data reduction are discussed 
in Sect. [2] We combine the H2O upper limits with published CO 
7=1-0 data to derive upper limits to the H2O abundance rela- 
tive to H2, x(o-H20), in Sect. [3] We summarize and discuss the 
possibilities for future H2O detections with the Herschel Space 
Observatory in Sect. [4] 



2. Observations and Analysis 

We observed the central regions of six nearby starburst galaxies 
(NGC253, IC342, M82, NGC4258, CenA, and M51) with Odin 
between 2002 November 22 and 2005 June 17 with on-source 
integration times ranging from 17 to 25 hours. Odin's 1.1m tele- 
scope has a beam diameter of 2.1' at the 556.936 GHz frequency 
of the 1 10 — loi transition of ortho-fyO and a main beam effi- 
ciency of 0.9 (Frisk et al. 120031 ). Details of the observations for 
each galaxy are given in TableQ] The observations were made in 
position switching mode; the off positions were located 850" to 
the south-east for NGC253, IC342, M82, and NGC4258, 280" 
to the south-east for CenA, and 900" to the east for M5 1 . The 
backend for all observations was an acousto-optic spectrometer 
with a spectral resolution of 1 MHz (0.5 km s~') and a total 
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Fig. 1. Spectra of nearby starburst galaxies tuned to the 1 iq - Iqi transition of ortho-WtO at 556.936 GHz. The velocity scale is 
Vlsr- No broad line (consistent with the observed CO 7=1-0 emission from these sources) is seen in any of the spectra, (a) NGC253 
(Vlsr = 249 km s~'). The spectrum was obtained toward coordinates 00:47:35.2 -25:17:20 (J2000); this position is 28" from the 
starburst nucleus, (b) IC342 (Vlsr = 32 km s~'). (c) M82 {Vlsr = 250 km s~'). Due to drifts in Odin's pointing model with time, 
this spectrum was obtained at an offset of ~ 55" from the center of M82, which is at 09:55:54.0 +69:40:57 (J2000). (d) NGC4258 
(V L sr = 472 km s" 1 ). (e) CenA (Vlsr = 550 km s" 1 ). (f) M51 (Vlsr = 460 km s" 1 ). 



bandwidth of 1.1 GHz (~ 600 km s~'). Odin's pointing is esti- 
mated to be accurate to 15". However, for M82, maps of Jupiter 
made after the data were obtained suggest that the pointing for 
this galaxy was in error by ~ 55". We take this pointing error 
into account in our analysis below. 

The data reduction followed the methods described in 
Wilson et al. (2003) to identify and remove individual spectra 
with anomalous system temperatures or contamination by the 
Earth's atmosphere. The individual spectra were averaged and 



the resulting spectrum was binned to a spectral resolution of 2.7 
km s~ l for the final analysis. The summed spectra were fit with 
a 0-2 order polynomial baseline excluding the expected velocity 
range of the H2O emission; typically a velocity range of 200-400 
km s -1 was excluded from the baseline fit. The resulting spectra 
are shown in Figure Q] 

Upper limits (3cr) to the integrated intensity of the H2O line 
were calculated using the width of the CO 7=1-0 line observed 
over the same area as the Odin beam. The upper limits to the 
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Table 1. H2O Observations and Abundance Limits 



Source 


Coordinates 


1 sys 


fill! 


rms fl 






U 


A/1 C 


X{0-ti2U)_ 


CO Reference 




(J2000) 


(K) 


(hr) 


(mK) 


(km s" 1 ) 


(Kkms" 1 ) 


(Mpc) 


(10* Mo) 






NGC253 


00:47:35.2 -25:17:20 


4100 


18.5 


13 


450 


1.6 


2.5 


0.60 


< 2.0 x 10-* 


Sorai et al. 2000 


IC342 


03:46:49.6 +68:05:45 


3700 


25.0 


10 


130 


0.56 


3.9 


0.39 


< 2.6 x 10-' 


Heifer et al.l2003l 


M82 


09:55:54.0 +69:40:57 


3400 


20.1 


10 


400 


1.9 


3.9 


1.56 


< 1.7 x 10- 9 


Walter etal.|2002l 


N4258 


12:18:57.5 +47:18:14 


4300 


17.4 


13 


450 


1.4 


8.1 


0.81 


< 1.3 x 10- 8 


Heifer et al. [20031 


CenA 


13:25:27.6 -43:01:08 


4500 


17.2 


10 


600 


1.2 


4.0 


0.29 


< 7.8 x 10-» 


Eckartetal. 179901 


M51 


13:29:53.1 +47:11:48 


3900 


23.7 


8 


200 


0.56 


7.7 


1.64 


< 2.4 x 10~ 9 


Heifer et al. [20031 



" Rms noise at a resolution of 2.7 km s 1 

* Three sigma upper limit integrated over AV(CO); values for M82 and NGC253 include scaling by 1.61 and 1.15, respectively, to account 
for pointing errors; see text 
Assumes X co = 2 x 10 20 H, cnT 2 (K km s^ 1 ) -1 

d Three sigma upper limit assuming n Hl = 10 6 cirr 3 and T K = 40 K. Depending on the average density of the molecular gas on large scales 
and the mass fraction in dense cores, these abundance limits could be larger by a factor of ~ 10 or more (see text). 



H2O integrated intensity and the CO line widths are given in 
Table [JJ CO line widths (full width at zero intensity) appropri- 
ate to the Odin beam were calculated for three galaxies from the 
data in the BIMA SONG survey (Heifer et al. 120031 . For C enA 
and M82, the emission in the maps of Eckart et al. (1990} and 
Walter et al. ([2002 ), respectively, lies almost entirely within the 
Odin beam and so we used all the emission to estimate the line 
width. For NGC 253, the full width of the CO line is de nned 
by the spectrum at the starburst nucleus (Sorai et al. 2000). For 
M82, the integrated intensity is corrected for the pointing error 
assuming a compact source and an Odin primary beam of 2.1' 
(full width at half maximum). For NGC253, the observed coor- 
dinates are 28" from the peak emission in the galaxy, and so we 
also correct the integrated intensity upper limit for this offset. 

3. Upper Limits to the H 2 0/H 2 Abundance Ratio 

We can use the upper limits to the FbO integrated intensity to 
calculate upper limits for the H2O abundance within the Odin 
beam using the formula given by Snell et al. (2000b). We adopt 
an H2 volume density of 1 x 10 6 cirr 3 and a temperature of 40 
K to calculate the H2O column density upper limits. We adopt a 
relatively high H2 volume density in these calculations because 
the H2O emission most likely originates in the densest parts of 
the interstellar medium. 

We use published CO 7=1-0 data to calculate the mass 
of molecular hydrogen within the Odin beam. For IC342, 
NGC4258, and M51, we integrated the CO flux within the Odin 
beam from the publicly available BIMA SONG data (Heifer et 
al. |2003l . We adopt a CO-to-H 2 con versio n factor X co = 2x 10 20 
H2 cirr 2 (Kkms -1 ) -1 (Strong et al. 1 19881 1 and use the formula in 
Wilson & Scoville (1990) to calculate the molecular gas mass. 
For CenA, we take the molecular gas mass from Eckart et al. 
( 1990) and rescale it to our adopted conversion factor and dis- 
tance. For M82, we take the molecular gas mass from Walter et 
al. (120021 but adjust the mass in the streamers to our adopted 
conversion factor. For NGC253, we adopt the mass in the bar 
plus the "nuclear bar" from Sakai et al. 2000 and scale it to our 
adopted distance and conversion factor. The molecular hydrogen 
gas masses and distances are given in Table [TJ along with the re- 
sulting upper limits for the ortho-WtO abundance, x(o-H20). 

Our limits on the ortho-WzQ abundance range from 2 x 10 -9 
to 1 x 10~ 8 . Note that these estimated abundance limits are in- 
versely proportional to our assumed value of 10 6 cm -3 for the H2 
volume density. The tightest constraints are for the four nearest 



and/or most gas-rich galaxies: NGC 253, IC342, M82, and M5 1 . 
For NGC4258, its larger distance is the primary reason why the 
abundance limit is larger than the other three galaxies, since its 
molecular hydrogen mass translates into a smaller average col- 
umn density in the Odin beam. For CenA, the smaller molecular 
hydrogen column density is also the primary reason for its larger 
abundance limit. 

Goicoechea et al. (2005) have detected H2O in absorption 
towards NGC 253 using the Infrared Space Observatory (ISO). 
They derive a lower limit to the H2O column density of > 1 x 10 15 
cm" 2 . This limit is significantly larger than the 3cr upper limit to 
the H2O column density derived from our Odin data, which is 
< 4 x 10 13 cm -2 . This apparent discrepancy can be understood 
by examining the regions probed by the two instruments: the 
Odin upper limit is an average over a 2.1' region, while the ISO 
absorption spectrum probes the much smaller region of the cen- 
tral continuum source. The central source has a size of ~ 12" 
at mid-infrared wavelengths (Forster Schreiber et al 2003) while 
examining archival 850 fim continuum images from the James 
Clerk Maxwell Telescope suggests the size of the central sub- 
millimeter source is no more than 30". These sizes correspond 
to filling factors of 0.01-0.06 in the larger Odin beam, which 
are sufficient to bring the two column density limits into agree- 
ment. Alternatively, assuming a density of 10 4 cm -3 (Bayet et 
al. 2004, Glisten et al. 2006) would also bring the two abundance 
estimates into good agreement (see below). 

These water abundance upper limits correspond to values av- 
eraged over kiloparsec scales that, in our own Galaxy, would in- 
clude many giant molecular clouds. If we were to picture the 
molecular interstellar medium in these starburst galaxies as be- 
ing similar in structure to normal spiral galaxies, then only a 
fraction of the molecular gas would be contained in cores that 
are dense enough to produce significant H2O emission. In this 
situation, the upper limit to the water abundance in the cores 
themselves would be larger than the values given in TableQ] per- 
haps by factors of 10-50 depending on the exact fraction of gas in 
dense corefl Alternatively, the molecular interstellar medium of 
starburst galaxies is likely to have considerably higher average 



1 For example, the W3 cloud contains a total of 10 5 M G (Lada et al. 
[1978), of which 1900 M Q (Tieftrunk et al. |1995t is contained within the 
single core studied in H2O emission by Wilson et al. (2003). Depending 
on the masses of the three other large cores in W3(OH) and W3 Main 
(Tieftrunk et al. |1995|[T998l , this cloud likely contains 2-8% of its mass 
with densities high enough to excite water emission. 
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densities than those of normal galaxies. If water emission arises 
from an extended medium with a uniform density of 10 4 - 10 5 
cm -3 (Bayet et al. 120041 Giisten et al. 2006 ), the upper limits on 
large scales given in TableQ]would increase by factors of 10-100. 

It is interesting to compare these upper limits to the H2O 
abundance in starburst galaxies with the H2O abundance esti- 
mated for Galactic star forming regions. Snell et al. (2000a|) es- 
timate H2O abundances ranging from 6 x 10" 10 to 1 x 10" 8 for 
eight giant molecular cloud cores. (They also obtain upper limits 
for the H2O abundance in dark clouds, with the most stringent 
upper limit of 7 x 10~ 8 being obtained for TMC-1.) Wilson et 
al. (|2003 ) obtain an H 2 abundance for the W3 IRS5 region of 
2 x 10" 9 . Olofsson et al. d2003l used Odin to detect H^O at two 
positions in Orion and derive an H2O abundance of ~ 10~ 5 to- 
ward the outflow associated with Orion BN-KL and (1 - 8)x 10~ 8 
in the ambient cloud 2' to the south. More recently, Wirstrom et 
al. (|2P06) have used C 18 7=5-4 and H2O observations from 
Odin to measure a water abundance of > 8 x 10~ 8 in the Orion 
photon-dominated region. Plume et al. (2004) have used absorp- 
tion line measurements toward W49A to obtain H2O abundances 
of 8 x 10~ 8 to 4 x 10~ 7 toward three foreground clouds. Bergin 
et al. ( 120031) found much higher abundances of > 10~ 6 toward 
three outflows in the NGC 1333 molecular cloud core. Snell et 
al. (120051) have detected H2O emission in the supernova remnant 
IC443 with abundances (assuming [CO]/[H2] = 10~ 4 ) between 
2 X 10~ 8 and 3 x 10~ 7 . They suggest that both photodissociation 
of H2O in post-shock gas and freeze-out of H2O in the ambient 
gas are needed to explain the weakness of the H2O lines. 

The above discussion shows that the Odin upper limits for 
the H2O abundance averaged over kiloparsec scales are compa- 
rable to the lowest abundances obtained in the cores of Galactic 
giant molecular clouds. Naively, one might expect the interstel- 
lar medium in intense starburst galaxies like M82 to have simi- 
lar physical conditions to high-mass star forming cores but with 
those conditions now present on kiloparsec scales. The H2O up- 
per limits obtained here suggest that this naive picture is not cor- 
rect, and that even in intense starburst galaxies there is a sub- 
stantial mass of molecular gas that is not in the warm dense 
conditions appropriate to produce H2O emission. This conclu- 
sion is consistent with recent estimates of the average gas den- 
sity of 10 4 - 10 5 crrT 3 in some nearby starburst galaxies (Bayet 
et al. 120041 Giisten et al. 12006b . It is also possible that the in- 
tense ultraviolet radiation fields produced in starburst galaxies 
photo-dissociate more of the H2O molecules than is the case 
for individual high-mass cores like W3. Melnick et al. (2005 ) 
suggest that H2O in the gas phase is restricted to a thin layer 
of gas in molecular cloud cores, deep enough that it is shielded 
from the ambient ultraviolet radiation field but not so deep that 
the molecules freeze out onto dust grains. The stronger ultravi- 
olet radiation fields in starburst galaxies might conceivably al- 
ter this delicate balance and make H2O an even rarer molecule 
in the gas phase than it is in our own Galaxy. Finally, dramatic 
H2O absorption features are seen in observations toward Sgr A 
(Sandqvist et al. 120031 and Sgr B2 (Neufeld et al. 120031 in the 
Galactic Center region. These observations may imply that self- 
absorption of the H2O line could make it very difficult to detect 
emission in regions such as starburst nuclei with warm dense gas 
containing multiple outflows and shocks from supernova rem- 
nants. All these effects may conspire to make the emission lines 
of H2O very difficult to detect even in the brightest and closest 
starburst galaxies. 



4. Conclusions and Future Prospects 

We have obtained long integrations with the Odin satellite in an 
attempt to detect the emission line of orf/10-water in six nearby 
starburst galaxies. None of the galaxies were detected, with up- 
per limits to the H2O integrated intensity ranging from 0.6 to 1 .9 
K km s . We have combined these upper limits with published 
CO data to derive upper limits to the H2O abundance in each 
galaxy. Our 3<x upper limits to [FhOJ/tFb] range from 2 x 10~ 9 
to 1 X 10~ 8 . The most sensitive limits are comparable to the mea- 
sured abundance of H2O in the Galactic star forming region W3 
(Wilson et al. 2003). However, if only 10% of the molecular gas 
is in very dense cores, then the water abundance limits in the 
cores would be larger by a factor of 10 i.e. 2 x 10~ 8 to 1 x 10" 7 . 
Similarly, if the average density in these starburst galaxies is 
10 4 - 10 5 crrT 3 , the abundance limits would be larger by fac- 
tors of 10-100. 

These observations can provide useful guidance for planning 
additional searches for H2O emission with future space mis- 
sions. The most sensitive 3cr upper limit to the H2O emission 
obtained with Odin is 0.5 K km s -1 (for IC342 and M51), which 
corresponds to 1600 Jy km s _1 or 3 x 10~ 17 W irT 2 . The cur- 
rently quoted sensitivity for the HIFI 460-680 GHz heterodyne 
spectrometetQ on the Herschel Space Observatory corresponds 
to a 3cr upper limit of 1.4 x 10~ 18 W irT 2 for a spectral line of 
width 130 km s for one hour of dual beam switched obser- 
vations, or about twenty times more sensitive than our existing 
Odin data. Thus, to obtain H2O spectra of these galaxies which 
are ten times more sensitive than the Odin spectra presented here 
will require integration times with Herschel of only 5-20 minutes 
per source. 

Alternatively, we can estimate what integration time will be 
needed to reach an H2O abundance limit of 2 x 10~ 10 , e.g., a 
factor of 10 lower than the abundance of water in W3 (Wilson 
et al. 120031 . For NGC253 and M82, this limit is reached in 3-6 
minutes, while IC342 and M5 1 would require ~ 30 minutes of 
integration. However, CenA and NGC4258 require 3.5 and 5.1 
hours, respectively, to reach this abundance limit. These longer 
integration times are primarily due to the lower surface density 
of molecular gas in these two galaxies. This analysis suggests 
that significant integration times may be required to detect H2O 
in emission even with Herschel, and that the best targets for 
the initial searches will be the most gas-rich, nearby starburst 
galaxies. If deep searches are not successful in detecting H2O 
in emission, it may have interesting implications for the astro- 
chemistry of starburst galaxies. However, the recent ISO results 
on H2O absorption lines in starburst galaxies (Goicoechea et al. 
2005, Gonzalez-Alfonso et al. 2004) suggest that absorption line 
searches with Herschel may be very fruitful, especially given its 
much improved angular resolution compared to Odin. 
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